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Abstract 

We consider interference patterns produced by coherent arrays of Bose-Einstein 
condensates during their one-dimensional expansion. Several characteristic pattern 
structures are distinguished depending on value of the evolution time. Transformation 
of Talbot "collapse-revival" behavior to Fraunhofer interference fringes is studied in 
detail. 

1 Introduction 

The interference measurements ^ Ej on two expanding Bose-Einstein condensates (BECs) 
have created new important field of research where the density profile of gas, imaged after 
releasing from the trap, provides important information about the phase of the ground- 
state wave function. Expansion of coherent arrays of BECs provides new opportunities 
to test the phase properties of the system For example, Fraunhofer interference 

patterns observed in [3j demonstrate strong coherence of BECs confined in separate traps, 
and experiment [U] shows that this coherence can be manipulated by means of collapses and 
revivals of wave functions due to nonlinear interaction of BECs in tightly confined states 
formed by three-dimensional periodic trapping potential. 

It is well known (see, e.g., [?]) that mentioned above "collapse- revival" behavior of 
quantum-mechanical wave functions jHl El is a temporal counterpart of optical Talbot ef- 
fect fOl ^] in which interference pattern behind the grating restores at distances multiple 
of the so-called Talbot distance /\ {d is the slit spacing in the grating and A is the wave- 
length of light). Similar Talbot effect was also observed in atom optics jil2| I13j. Analogy 
between spatial Talbot effect and temporal collapse-revival behavior of wave functions sug- 
gests that such collapses-revivals should exist in interference of matter waves emitted from 
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arrays of BECs provided evolution time is small enough, and indeed such effect was observed 
in Plj. In this connection it is natural to ask how this short-time Talbot behavior evolves 
for finite array of condensates to long-time Fraunhofer behavior observed in This paper 
is devoted to consideration of this problem. 

In Section 2 we present general formulas for the wave function produced by a linear 
array of BECs. We confine ourselves with one-dimensional theory under supposition that 
condensate remains confined in radial direction after turning off a periodic optical potential 
and evolution takes place only along the axial direction of the BECs array. This formulas 
permit us to distinguish characteristic stages of evolution — short-time Talbot stage with 
revivals of the wave function in the central part of the array, intermediate time stage when 
Fraunhofer fringes already formed with Fresnel diffraction pattern inside each of them, and 
long-time Fraunhofer stage with standard density distribution along fringes. These stages of 
evolution of the wave function are studied in detail in Section 3 (Talbot stage) and Section 
4 (transition to Fraunhofer stage). The last Section 5 is devoted to conclusions. 



2 General formulas 



After switching off the periodic optical potential the condensate density decreases and under 
condition that the initial size ~ cr of each BEC is much less than the spacing d between sites, 
the interatomic interaction can be neglected during most time of the evolution and, hence, 
the wave function obeys the linear Schrodinger equation 



ihipt = -—ipxx- 
2m 

If the initial state is given by ip{x, 0) = ipo{x), then after time t it evolves into 



ip{x,t) = / G{x — x' ,t)4'o{x')dx' 

J — oo 

where G{x — x', t) is well-known Green function of Eq. 

G{x-x',t) = 



(1) 



(2) 
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see, e.g. [T^): 

/\2" 



2ht 
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To simplify calculations, we suppose that the initial wave function of BEC in the site of the 
array with the coordinate kd can be approximated by a Gaussian function and, hence, we 
represent the initial state of BEC as 



[x — kdy 
2^2 



(4) 



where Nk = l^fcP is equal to number of atoms in kth condensate (we suppose that a <^ d) 
and (pk is its phase. Then Eq. 1^ yields the solution 



il){x,t) 
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exp 



{x - kdy^ 



2a'^{l+iht/ma^)_ 



(5) 
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This formula should be specified in accordance with the problem under consideration. In the 
case of large number of condensates in the array confined in axial direction by a parabolic 
potential, the Thomas-Fermi approximation can be used for calculation of number of atoms 
in kth condensate which gives jl]: 

Nk = ^l = T^^-T:r) ' (6) 



16kM V 
where 



2hjuj f 15 a d 



_ 1/5 



= is the total number of atoms, uj = {uJxUj'j_Y^^ is the geometric mean of the 

magnetic trap frequencies, a/jo = ^JfLfm^ is the corresponding oscillator length, and a > 
is the s-wave scattering length. 

In the experiment |3j there was ku — 10^ ^ 1, and this large parameter suggests that 
there are different stages of evolution. 

For 

md^ 

t « — (8) 

each condensate evolves independently of each other and there is no their interference effects. 
For 

md'^ , md'^ 

t ^ < A;m ■ ^ 9 

n h 

we have interference between condensates, but in the central part of the array the infiuence 
of its finite size is negligibly small and local interference pattern can be approximated by 
that of an infinite periodic lattice of condensates which leads to temporal Talbot effect. 
For 

kM-^<^t<t:klj-—- 10 
n h 

the Fraunhofer fringes begin to form. Indeed, their positions are given by (see, e.g., |3]) 

x^{t) = ±n^t, n = 0, 1, 2,..., (11) 
dm 

and if t satisfies the condition (fTUI) . then distances between neighboring fringes ~ 27cht/md 
are much greater than the size of each fringe ~ 2kMd (see below). At the same time, the 
interference pattern inside each fringe is formed by only some part of the array and hence 
we get Fresnel diffraction pattern along the fringe. 
At last, for 

t » kl ■ -J- (12) 

we arrive at usual Fraunhofer diffraction when the whole array contributes into interference 
pattern inside each fringe. 

To illustrate these stages of evolution of the wave function, we have shown in Fig. 1 the 
distributions of density p = l?/'^ calculated from formulas ©-((Zj) with 0^ = (coherent con- 
densates). We see that for t <^ kMijnd"^ /h) the density distribution reproduces periodically 
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Figure 1: Evolution of the density profile of array of condensates with time calculated 
according to Eqs. ©-((Tj) with ci = 8, ci = 1 (in dimensionless units) and /cm = 10. At 
t = O.Str, where is given by Eq. ((201); see complex interference pattern ("collapse" of 
wave function); at t = 0.5tr the central part coincides with that for t = but shifted to a 
half-period d/2; at t = tr the initial distribution is almost completely restored; at t = 5tr the 
side fringes start to form, and, finally, at t = 30tr we see Fraunhofer diffraction of matter 
waves from finite "grating" . 



in time with period tr ~ mdP /h (see exact formula ()20|1 below), the side fringes begin to form 
at t ~ kuijnd? /h), and for t S> kMimcf /h) there are peaks of density at the coordinates 
given by (fTT|) and profiles of fringes take Fraunhofer form for t ^ {kMd)'^m/h. The solution 
(0) permits us to investigate these stages of evolution analytically. 



3 Talbot revivals of wave function 

For time values in the region Q, we can approximate the array by infinite lattice of equidis- 
tant condensates so that for coherent condensates their wave function is given by 



ip{x,t) 



A 



oo 

^ exp 



[X 



kdf 



(13) 



4 



With the use of definition of ^3-function (see, e.g. fT^) 

00 

e3{z,T) = ^ exp[i7i{Tk^ + 2zk)] (14) 



fc=— 00 

the wave function p3p can be presented in the form 



TT 2ai ( i'nx^\ „ f x 1 



^,(.,*) = _W_exp 93K;,--1, (15) 



where 

2TTia'^ ( iht 



By means of transformation formula (see ^H] 



1 + Z-? • (16) 



^^3( V^) =v^exp('^)e3(^,r) (17) 



we transform (jlSp into 



,, , VSttct a ^ (X \ 



Then the periodicity property of 6'3-function, 6*3(^5 r ± 2) = ^^3(z,r), leads at once to the 
periodicity of the wave function, 

V^(x,t + t,) = 7/'(x,t), (19) 

with the period 

tr = —T- (20) 

nn 

If the array of BECs is realized in optical periodic potential with light wavelength A, then 
the spacing between neighboring lattice sites is equal to d = A/2 and the revival time can 
be expressed in the form 

= 4 ■ ^' 

where 

Er = ^ (22) 
"2m ^ ^ 

is the recoil energy (g = 2it/\). 

The evolution time tr/2 corresponds to the transformation of ^3-function 63{z,t + 1) = 
^4(^5 r) = ^3(-2 + 1/2, r), that is we obtain the wave function shifted to the distance d/2 with 
respect to its initial form: 

ilj{x,t + tr/2) =iIj{x + d/2, t). (23) 

Above calculation explains periodic restoration of initial wave function by means of trans- 
formation properties of ^-functions. To relate this approach with standard one (see, e.g., 0), 
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Figure 2: Evolution of density profiles for BEC arrays with zero relative phase. Two first 
revivals at t = ntr, n = 1,2, are clearly seen as well as "fractional revivals" at intermediate 
moments tr/S, tr/4:, 3tr/S, tr/2, etc. 



let us consider the problem from a different point of view. The linear Schrodinger equation 
with periodic initial condition 



^(x,0) 



A 



71 



1/4 



exp 



(x - kdy 
2^ 



(24) 



can be solved by the Fourier method which gives 




2„2 



27rV 



1 + 



iht 



2nx 

cos I 



(25) 



(The relationships between presentations (fT^ and of the same wave function ip{x,t) is 
expressed by the known identity for these two series; see, e.g., [HI.) We see that at t equal 
to multiple of Talbot time t,., t = ntr, all phase factors in ()25|) become equal to unity and 
(|25p reduces to the Fourier series for the initial periodic wave function (j24j) . This method of 
derivation of time-periodicity of the wave function shows that periodic restoration of the ini- 
tial state is not a specific feature of the initial state (j^^ built of Gaussian functions. Indeed, 
any periodic initial function can be expanded into Fourier series and harmonics cos{27Txk/d), 
fc = l,2,. 
unity at t 



evolve with time according to factors exp(— i^^/c ) which become equal to 



ntr- Thus, any periodic initial wave function completely restores periodically 
its form. The described above picture of periodic in time changes of the interference pattern 
is shown in Fig. 2 where even for relatively small number of condensates first several revivals 
are clearly seen. 
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The above theory can be generahzed on non-zero phases in the initial state and hence in 
the solution (j21). For example, in the case of alternating phases of condensates, 

gi<^fc = (-1)^, (26) 
the wave function can be expressed in terms of 6'4-function [T^ . 



l2ai ( i-nx^ \ ^ I x 



^.fe*) = _y_exp --^ 1. (27) 



or, with the use of the transformation formula [T^ . 



e, ( i,-!") = v/=z7exp ) 02(^,r), (28) 



r r 



in the form 

Then the property 92{z,t + 1) = exp(7ri/4)6'2(z, r) leads to restoration of the initial state 
(up to inessential constant phase factor) after revival time 

md"^ 1 2TTh , . 

tr = — r = T^- 30 

Let us estimate an order of magnitude of the revival time for arrays of BECs. In the 
case [U] of '^^Rb BECs array loaded into optical potential with light wavelength A = 838 nm 
formula (PT|) gives tr — 75 /xs. This is about one order of magnitude less than the revival time, 
caused by nonlinear interaction, of single condensate in the experiment [Hj. In this experiment 
absorption images were taken after a time-of-fiight period of 16 ms which is much greater 
(with factor ~ 200) than our estimate of t^. For number of sites in 3D lattice ~ 10^ we have 
kM ~ 10 and, hence, the observed interference patterns correspond to the Fraunhofer limit 
()12|). In this case the difference of interference patterns was caused by difference in initial 
states of condensates at different "hold times" of evolution of each condensate in strongly 
confined states formed by 3D periodic trapping potential. 

In the experiment |3] the revival time is tr — 69 /is and a typical image was taken at 
t = 29 ms, that is for k^j ~ 100 again in the Fraunhofer limit (in accordance with the theory 
developed in this paper). 



4 Transition to Fraunhofer interference 

Now we shall turn to the regions ()10p and ()12j). Effects of Fresnel diffraction can be noticed 
in Fig. 1 for t = lOtr- However, they are not expressed clearly enough because of smooth dis- 
tribution (jni) of density in the array used in our calculations. Therefore it is more instructive 
to consider finite array with equal amplitudes Ak = 1 of wave functions in each condensate 
and take (pk = ■ k, that is with equal differences A0 of phases between neighboring 
condensates. Then Eq. (jSj) with t ^ ma'^/h reduces to 

1 I ma imx^ _iB££^ sr^ 
^(^'^) = ;^A/-^e2"e ^ ^ exp 

k=-kM 



mdx 

ht 



A(f) k 



irad 
2ht 



(31) 
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where we have taken into account only leading real and imaginary contributions in the series 
expansion of the exponential in powers of ma^/ht. The sum here has maximal amplitude 
when all terms are in phase in linear in k approximation. This condition defines coordinates 
Xn of the centers of fringes, 



2'Kh 
md 



nH -]t, n = 0, ±1,±2, ... 

2tc J 



(32) 



To consider profiles of fringes, we introduce the coordinate 6 which is reckoned from the 
center of the fringe: 

X = Xn + S, (33) 

so that dependence on S is determined mainly by the factor 



$((5,t) = J2 



k=—kf. 



imd6 , imd"^ , , 
-k + —^k^ 



ht 



2ht 



(34) 



If t satisfies the condition ()1U|). then both terms in the exponential have the same order 
of magnitude and, on one hand, the fringe width is of order of magnitude of the array 
length, 6 ~ 2kMd, and, on the other hand, it is much less than the distance between fringes. 
Therefore the coordinate x in the factor exp {—m'^a'^x'^ /2h'^t'^) can be replaced by Xn- Thus, 
the wave function in vicinity of the nth fringe is given by 



1pn{x,t) 



1 



ma 



1/4 V iht 



e 2fit ■ exp 



27rV2 / A0\' 



(35) 



where 6 = x — Xn- Now, for kM ^ 1 the sum in (jH^ can be approximated by integrals which 
are easily expressed in terms of Fresnel functions 18j: 



g 2ht 



C 



m 



2ht 



— {kMd + 6) +C \ —ikMd - 6) 



m 



2ht 



m 



Thus, distribution of density in the nth fringe is given by 



^(A:Mrf-5) 



(36) 



TTCT 



(i2 



exp 



X 



47rV 



n + 



A0\ 



2tx I 



m 



C{sl^^{kMd + 6) 



C 



^,ikMd-6) 



(37) 



+ 



^(W + 5) 



m 



^(W-5) 



The exponential factor determines the number of atoms in the nth fringe: 

N„, = const ■ exp P + 2^ J 



(38) 
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Figure 3: The central fringe profile for several values of the number of sites in the array. 
Time t corresponds to the region (fTUj) . The plots are calculated for d = 8, cr = 1 and (a) 
A;m = 20 at t = 40t^; (b) A;m = 40 at t = SOU; (c) fc^ = 80 at t = 160^. Formation of the 
Fresnel pattern is clearly seen. 

This formula reduces to Eq. (6) of Ref. |3] for A0 = 0. 

Dependence on 6 determines fine interference pattern inside fringes. It is expressed by 
the factor in curly brackets and demonstrates typical Fresnel form (see, e.g., [1^], section 
3.3, or [ini, section 8.7) of diffraction from a slit with width 2kMd equal to the whole array 
length. Accuracy of this analytical description depends on the number of sites in the array 
and increases with growth of /cm- In Fig. 3 it is shown how exact profile of density along 
the fringe changes with increase of /cm- Its transformation into Fresnel diffraction pattern is 
clearly seen. Small "ripples" are obviously caused by the discrete structure of the array. 

For larger values of time ()12|1 Fresnel structure evolves into usual form of density distri- 
bution in Fraunhofer diffraction from finite slit with width 2k Md. In this limit of very large 
t the quadratic in k tern in exponentials in Eq. (jMj) is much less than unity and can be 
omitted. Then simple integration gives ^{S,t) = sin (^^^^^) and hence distribution of 
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Figure 4: Evolution of density profile of central fringe on time. Values of the parameters 
are equal to d = 8, a = 1, ku = SO and (a) t = 500^; (b) t = lOOOt^; (c) t = 2000t^; (d) 
t = 16000tr-- Transformation of Fresnel profile shown in Fig. 3 (c) to standard Fraunhofer 
profile is clearly seen. 



density inside fringes is proportional to 



2 sin [mkMdd /ht) 



which is standard Fraunhofer diffraction distribution from finite slit (see, e.g., ^H]; section 
8.5). The described here evolution of profile is illustrated in Fig. 4. The total intensity of 
nth fringe is still determined, of course, by Eq. 



5 Conclusion 

We have presented in this paper analysis of interference of matter waves during one-dimensional 
expansion of finite arrays of condensates. It shows that the interference pattern exhibits 
quite complicated evolution with time from Talbot "collapses and revivals" of wave function 
through intermediate region of Fraunhofer fringes with Fresnel patterns inside them, and, 
eventually, to standard Fraunhofer diffraction from finite grating. One may suppose that 
technique of density imaging will permit one to study experimentally all these stages. 
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